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ABSTRACT 

We have compiled a complete, significance limited extragalactic sample based on ~ 
25, 000 deg 2 to a limiting flux of 3 x 1CT 11 ergs cm -2 s _1 (~ 7, 000 deg 2 to a flux limit of 
10 -11 ergs cm -2 s -1 ) in the 20 - 40 keV band with INTEGRAL . We have constructed a detailed 
exposure map to compensate for effects of non-uniform exposure. The flux-number relation is 
best described by a power-law with a slope of a = 1.66 ± 0.11. The integration of the cumu- 
lative flux per unit area leads to /20-40keV = 2.6 x 10 -10 ergs cm” 2 s -1 sr -1 , which is about 
1% of the known 20 - 40keV X-ray background. We present the first luminosity function of 
AGN in the 20-40 keV energy range, based on 68 extragalactic objects detected by the imager 
IBIS/ISGRI on-board INTEGRAL . The luminosity function shows a smoothly connected two 
power-law form, with an index of 71 = 0.9 below, and 72 = 2.2 above the turn-over luminos- 
ity of L* — 4.6 x 10 43 ergss“ 1 . The emissivity of all INTEGRAL AGNs per unit volume is 
W20— 40 keV (> 10 41 ergss” 1 ) = 2.8 x 10 38 ergss _1 h 3 0 Mpc~ 3 . These results are consistent with 
those derived in the 2 - 20keV energy band and do not show a significant contribution by 
Compton-thick objects. Because the sample used in this study is truly local (z = 0.022), only 
limited conclusions can be drawn for the evolution of AGNs in this energy band. But the objects 
explaining the peak in the cosmic X-ray background are likely to be either low luminosity AGN 
(Lx < 10 41 ergs s'* 1 ) or of other type, such as intermediate mass black holes, clusters, and star 
forming regions. 

Subject headings: galaxies: active — gamma rays: observations — X-rays: galaxies surveys — galaxies: 
Seyfert 
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1. Introduction 

The Galactic X-ray sky is dominated by accret- 
ing binary systems, while the extragalactic sky 
shows mainly active galactic nuclei (AGN) and 
clusters of galaxies. Studying the population of 
sources in this energy range has been a challenge 
ever since the first observations by rocket borne 
X-ray detectors (Giacconi et al. 1962). At soft 
X-rays (0.1 - 2.4 keV) the ROSAT All-Sky Survey 
(RASS; Voges et al. 1999) has revealed an extra- 
galactic population of mainly broad line AGNs, 
such as type Seyfert 1 and quasars. Between a 
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few keV and ~ 1 MeV, no all-sky survey using 
imaging instruments has been performed to date. 
In the 2-10 keV range surveys based on signif- 
icant fractions of the sky have been carried out 
with ASCA (e.g. Uedaet al. 2001), XMM-Newton 
(e.g. Hasinger 2004), and Chandra (e.g. Brandt 
et al. 2001) and have shown that the dominant 
extragalactic sources are more strongly absorbed 
than those within the RASS energy band, thus 
more Seyfert 2 type objects are detectable here. 
For a summary on the deep X-ray surveys below 
10 keV see Brandt & Hasinger (2005). At higher 
energies the data become more scarce. The Rossi 
X-ray Timing Explorer ( RXTE) sky survey in the 
3 — 20keV energy band revealed 100 AGNs, show- 
ing an ever higher fraction of absorbed sources of 
about 60% (Sazonov & Revnivtsev 2004). A sim- 
ilar study at energies above 20 keV had not been 
possible, because a large field of view in combi- 
nation with sufficient sensitivity are required to 
study the AGN population on a relevant fraction 
of the sky. The International Gamma-Ray Astro- 
physics Laboratory ( INTEGRAL ; Winkler et al. 
2003) offers an unprecedented > 20keV collect- 
ing area and state-of-the-art detector electronics 
and background rejection capabilities. Notably, 
the imager IBIS with an operating range from 
20 - 1000 keV and a fully-coded field of view of 
10° x 10° enables us now to study a large por- 
tion of the sky. A first catalog of AGNs showed a 
similar fraction of absorbed objects as the RXTE 
survey (Beckmann et al. 2006a). 

Related to the compilation of AGN surveys in 
the hard X-rays is the question of what sources 
form the cosmic X-ray background (CXB). While 
the CXB below 20 keV has been the focus of many 
studies, the most reliable measurement in the 10 
- 500 keV has been provided by the High Energy 
Astronomical Observatory (HE AO 1), launched in 
1977 (Marshall et al. 1980). The most precise 
measurement provided by the UCSD/MIT Hard 
X-ray and Gamma-Ray instrument (HE AO 1 A- 
4) shows that the CXB peaks at an energy of 
about 30keV (Gruber et al. 1999). The isotropic 
nature of the X-ray background points to an ex- 
tragalactic origin, and as the brightest persistent 
sources are AGNs, it was suggested earlyon that 
those objects are the main source of the CXB (e.g. 
Zdziarski 1996). In the soft X-rays this concept 
has been proven to be correct through the obser- 


vations of the ROSAT deep X-ray surveys, which 
showed that 90% of the 0.5 - 2.0 keV CXB can be 
resolved into AGNs (Miyaji, Hasinger & Schmidt 
2000). At higher energies (2 - lOkeV), ASCA 
and Chandra surveys measured the X-ray lumi- 
nosity function (XLF). These studies show that 
in this energy range the CXB can be explained 
by AGNs, but with a higher fraction of absorbed 
(Nh > 10 22 cm~ 2 ) objects than in the soft X-rays 
(e.g. Ueda et al. 2003). A study based on the 
RXTE survey by Sazonov & Revnivtsev (2004) de- 
rived the hard X-ray luminosity function of AGNs. 
They showed that the summed output of AGNs in 
this energy range is too small to explain the CXB, 
and suggested that a comparable X-ray flux may 
be produced together by lower luminosity AGNs, 
non-active galaxies and clusters of galaxies. 

With the on-going observations of the sky by 
INTEGRAL , a sufficient amount of data is now 
available to derive the AGN hard X-ray lumi- 
nosity function. In this paper we present anal- 
ysis of recent observations performed by the IN- 
TEGRAL satellite, and compare the results with 
previous studies. In Section 2 we describe the 
AGN sample and in Section 3 the methods to de- 
rive the number-flux distribution of INTEGRAL 
AGNs are presented together with the analysis of 
their distribution. Section 4 shows the local lumi- 
nosity function of AGNs as derived from our data, 
followed by a discussion of the results in Section 
5. 

2. The INTEGRAL AGN Sample 

Observations in the X-ray to soft gamma- 
ray domain have been performed by the in- 
struments on-board the INTEGRAL satellite 
(Winkler et al. 2003). This mission offers the 
unique potential to perform simultaneous observa- 
tions over the 2 - 8000 keV energy region. This is 
achieved by the X-ray monitor (2-30 keV) JEM- 
X (Lund et al. 2003), the soft gamma-ray im- 
ager (20-1000 keV) ISGRI (Lebrun et al. 2003), 
and the spectrograph SPI (Vedrenne et al. 2003), 
which operates in the 20 - 8000 keV region. Each 
of these instruments employs the coded-aperture 
technique (Caroli et al. 1987). In addition to 
these data an optical monitor (OMC, Mas-Hesse 
et al. 2003) provides photometric measurements 
in the V band. 


2 


Data used for the analysis presented here were 
all in the public domain by the end of March 2005, 
This includes data from orbit revolutions 19 - 137 
and revolutions 142 - 149. Data before revolution 

19 have been excluded as the instruments settings 
changed frequently and therefore the data from 
this period are not suitable to be included into a 
homogeneous survey. 

The list of sources was derived from the analysis 
as described in Beckmann et al. (2006a). Addi- 
tional observations performed later led to further 
source detections within the survey area. We ex- 
tracted spectra at those positions from the data 
following the same procedure. It is understood 
that most of those objects did not result in a sig- 
nificant detection in the data set used here, but it 
ensures completeness of the sample at a moderate 
significance limit. 

The analysis was performed using the Offline 
Science Analysis (OS A) software version 5.0 dis- 
tributed by the ISDC (Courvoisier et al. 2003a). 
We applied the same method for IBIS/ISGRI and 
SPI analysis as described in Beckmann et al. 
(2006a). The analysis of the INTEGRAL /IBIS 
data is based on a cross-correlation procedure be- 
tween the recorded image on the detector plane 
and a decoding array derived from the mask pat- 
tern (Goldwurm et al. 2003). The ISGRI spectra 
have been extracted from the count rate and vari- 
ance mosaic images at the position of the source, 
which in all cases corresponds to the brightest 
pixel in the 20 — 40keV band. 

The significances listed in Tab. 1 have been de- 
rived by using the OSA software, and refer to the 
count rate and count rate error for ISGRI in the 

20 - 40 keV energy band. The luminosities listed 
are the observed luminosities in this energy band. 
The absorption listed is the intrinsic absorption 
in units of 10 22 cm -2 as measured in soft X-rays 
by various missions as referenced. We also include 
the most important reference for the INTEGRAL 
data of the particular source in the last column of 
Table 1. The extracted images and source results 
are available in electronic form 1 . 

In order to provide a complete list of AGNs 
detected by INTEGRAL , we included also those 
sources which are not covered by the data used 
for our study. Those sources are marked in Tab. 1 

1 http://heasarc.gsfc.nasa.gov/docs/integral/inthp_archive.htmI 


and are not used in our analysis. 

3. Number-Flux Distribution of INTE- 
GRAL AGNs 

In order to compute the AGN number-flux re- 
lation counts of AGNs it is necessary to have a 
complete and unbiased sample, and to know the 
characteristics of the survey used to extract the 
data. Because of the in-homogeneous character of 
the survey, we had to apply a significance limit 
rather than a flux limit to define a complete sam- 
ple. The task is to find a significance limit which 
ensures that all objects above the given limit have 
been found. To test for completeness, the V e /V a - 
test has been applied. 

The V e / P a -test is a simple method developed 
by Avni & Bahcall (1980) based on the V/V max 
test of Schmidt (1968). V e stands for the vol- 
ume, which is enclosed by the object, and V a is 
the accessible volume, in which the object could 
have been found (e.g. due to a flux limit of a 
survey). Avni & Bahcall showed that different 
survey areas with different flux limits in various 
energy bands can be combined by the V e /V a -test. 
In the case of no evolution (V e /V a ) = 0.5 is ex- 
pected. This evolutionary test is applicable only 
to samples with a well-defined significance limit 
down to which all objects have been found. It can 
therefore also be used to test the completeness of 
a sample. We performed a series of V e /V a -tests 
to the INTEGRAL AGN sample, assuming com- 
pleteness limits in the range of la up to 6a IS- 
GRI 20 - 40 keV significance. For a significance 
limit below the true completeness limit of the sam- 
ple one expects the V e /V a -tests to derive a value 
(V e /V a ) < (V e /V a ) true , where ( V e /V a ) true is the 
true test result for a complete sample. Above the 
completeness limit the (V e /V a ) values should be 
distributed around (V e /V a ) frue within the statis- 
tical uncertainties. An example of this method to 
determine the significance limit can be found in 
Beckmann et al. (2006a). 

It appears that the sample becomes complete 
at a significance cutoff of approximately 5a, which 
includes 38 AGNs. The average value is (V e /V a ) = 
0.45 ± 0.05. This is consistent with the expected 
value of 0.5 which reflects no evolution and an even 
distribution in the local universe. 

A pure number flux distribution (i.e. log iV >5 
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versus log 5) for the sample presented here would 
not give meaningful results, because of the differ- 
ent exposure times across the survey, and there- 
fore the varying sensitivity within the survey. An 
uncorrected number flux distribution for INTE- 
GRAL AGNs has been shown in Beckmann et 
al. (2006b). To correct for the different expo- 
sure times it is necessary to count the number 
of AGNs per unit sky area. Thus the number 
of AGNs above a given flux have to be counted 
and divided by the sky area in which they are de- 
tectable throughout the survey. We therefore first 
determined the exposure time in 64, 620 sky ele- 
ments of — 0.63 deg 2 size within our survey. In 
each sky bin, the exposure is the sum of each indi- 
vidual exposure multiplied by the fraction of the 
coded field of view in this particular direction. The 
dead time and the good time intervals (GTI) are 
not taken into account but the dead time is fairly 
constant (around 20%) and GTI gaps are very rare 
in IBIS/ISGRI data. Figure 1 shows the exposure 
map in Galactic coordinates for this survey. We 
excluded those fields with an exposure time less 
than 2 ks, resulting in 47,868 sky elements with 
a total coverage of 9.89 sr. The flux limit for a 
given significance limit should be a function of the 
square root of the exposure time, if no systematic 
effects apply, but this assumption cannot be made 
here. The nature of coded-mask imaging leads to 
accumulated systematic effects at longer exposure 
times. In order to achieve a correlation between 
the exposure time and the flux limit, we therefore 
used an empirical approach. We correlated the ex- 
posure time of each object with an equivalent flux 
correspoding to 5cr significance based on the flux 
and significance of each sample object. By using 
only the AGNs of our survey we are additionally 
assured to consider effects based on the spectral 
slope (for AGNs T ^ 2). The correlation was then 
fit by a smooth polynomial (Fig. 2). This function 
was then used to estimate the flux limit of each 
survey field. It has to be noted that the individ- 
ual flux limit of each survey field is not important, 
but only the correct distribution of those flux lim- 
its. The total area in the survey for a given flux 
limit is shown in Figure 3. 

Based on the flux limits for all survey fields, we 
are now able to construct the number flux distri- 
bution for the INTEGRAL AGN, determining for 
each source flux the total area in which the source 


is detectable with a 5a detection significance in the 
20 — 40 keV energy band. The resulting correlation 
is shown in Fig. 4. 

3.1. The Slope of the Number-Flux Dis- 
tribution 

We applied a maximum-likelihood (ML) algo- 
rithm to our empirical number-flux distribution 
to obtain a power-law approximation of the form 
N(> S) = K • S~ a . Our approach was based on 
the formalism derived by Murdoch, Crawford and 
Jauncey (1973), also following the implementation 
of Piccinotti et al (1982) . The latter involved mod- 
ification of the basic ML incorporated to facilitate 
handling of individual source flux-measurement 
errors. The ML method also involves the applica- 
tion of a Kolmogorov-Smirnov (K-S) test to eval- 
uate the goodness of fit, as detailed in Murdoch, 
Crawford and Jauncey (1973). Once the slope is 
determined, a chi-square minimization is used to 
deterimine the amplitude K. 

For this analysis, we used a sub-sample of 38 
sources for which the statistical significance of our 
flux determinations was at a level of 5a or greater. 
The dimmest source among this sub-sample was 
fx = 5.6 x 10“ 12 ergscm" 2 s~ 1 , and the brightest 
was fx = 3.2 x 10" 10 ergscm~ 2 s~ A . We derived 
a ML probability distribution, which can be ap- 
proximated by a Gaussian, with our best fit pa- 
rameters of a = 1.66 ± 0.11. A normalization 
of K — 0.44 sr -1 (10" 10 ergs cm -2 s" 1 )* was then 
obtained by performing a least- squares fit, with 
the slope fixed to the ML value. 

4. The Local Luminosity Function of AGNs 
at 20 - 40 keV 

The complete sample of INTEGRAL AGNs 
with a detection significance > 5a also allows us 
to derive the density of these objects in the local 
Universe as a function of their luminosity. In or- 
der to derive the density of objects above a given 
luminosity, one has to determine for each source 
in a complete sample the space volume in which 
this source could have been found considering both 
the flux limit of each survey field and the apparent 
brightness of the object. We have again used the 
correlation between exposure time and flux limit 
as discussed in the previous section in order to as- 
sign a 5a flux limit to each survey field. Then the 
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maximum redshift z max at which an object with 
flux fx would have been detectable in each sky 
element was used to compute the total accessible 
volume 

N o. 

' = E 4^ WZmaxAfx)] (1) 

i= 1 

with N being the number of sky elements in which 
the object would have been detectable and f 2* the 
solid angle covered by sky element i. We applied a 
cosmology with H 0 — 70 km s'" 1 Mpc“ 1 ( hjo = 1), 
k = 0 (flat Universe), Q matter = 0.3, and Ao = 0.7, 
although a A 0 = 0 and qo = 0.5 cosmology does 
not change the results significantly because of the 
low redshifts in our sample. Figure 5 shows the 
cumulative luminosity function for 38 INTEGRAL 
detected (> 5a) AGNs in the 20 - 40 keV energy 
band. Here the density <j> describes the number 

of objects per Mpc 3 above a given luminosity Lx : 
k 

(j> = ^ V a / with K being the number of objects 

i~ 1 

with luminosities > Lx- Blazars have been ex- 
cluded because their emission is not isotropic. The 
redshifts in the sample range from 2 = 0.001 to 
z = 0.129 with an average redshift of z = 0.022. 
Thus the luminosity function is truly a local one. 
Figure 6 shows the luminosity function in differen- 
tial form, indicating that the sample indeed does 
not show luminosity bins with large incomplete- 
ness compared to the rest of the sample. The 
errors are based on the number of objects con- 
tributing to each value. The differential XLF also 
shows, like the cumulative one, a turnover around 
Lx = (5 — 10) x 10 43 ergss -1 . 

Because our study is based solely on low red- 
shift objects, we are not able to constrain models 
involving evolution with redshift. Nevertheless we 
can compare the XLF presented here with model 
predictions from previous investigations. XLFs 
are often fit by a smoothly connected two power- 
law function of the form (Maccacaro et al. 1991) 

dfKL x ,z = Q) = \fLxV 1 (Lx V 2 1 ^ 

d log L x [\Lj \lJ . 

( 2 ) 

We fit this function using a least-squares 
method applying the Levenberg-Marquardt algo- 
rithm. The best fit values we obtained are A = 
(2.7±0.5) x 1CT 6 hi o Mpc" 3 , j x = 0.93±0.15, 72 = 


2.23±0.15, and L * = (4.6±2.0)xl0 43 /if 0 2 ergss" 1 . 

The errors have been estimated on the basis of the 
uncertainties of the density determination and of 
the flux measurement. These values are consis- 
tent with values derived from the 2 — lOkeV XLF 
of AGNs as shown by e.g. Ueda et al. (2003) 
and La Franca et al. (2005). For example the 
work by Ueda et al. (2003) reveals for a pure 
density evolution model the same values for A, 

71 , and 72 , but a higher L* = 1.29 x 10 44 ergs _1 . 

This ratio of 2.8 in the turnover luminosity can 
be easily explained by the different energy bands 
applied. A single power law with photon index 
of T = 1.9 in the range 2 - 40keV would lead to 
L[ 2 —io keV) / •^ / (20 — 40 keV) ~ 2.8, assuming no intrin- 
sic absorption. Thus it appears that the local lu- 
minosity function of AGNs in the 20 — 40keV band 
can simply be extrapolated from the 2 — lOkeV 
range. This has, of course, no implications for 
the XLF at higher redshifts. The values are also 
consistent with the luminosity function for AGNs 
in the 3 — 20 keV band as derived by Sazonov & 
Revnivtsev (2004) from the RXTE all-sky survey. 

Information about intrinsic absorption is avail- 
able for 34 of the 38 objects (89%) from soft X-ray 
observations. This enables us to derive the lumi- 
nosity function for absorbed (Nh > 10 22 cm“ 2 ) 
and unabsorbed sources, as shown in Figure 7. 

The absorbed sources have a higher density than 
the unabsorbed sources at low luminosities, while 
this trend is inverted at high luminosities. The lu- 
minosity where both AGN types have similar den- 
sities is about L( 2 0 - 4 OkeV) ~ 3 x 10 43 ergs“ 1 . This 
tendency is also evident when comparing the frac- 
tion of absorbed AGNs with the luminosity in the 
three luminosity bins depicted in Figure 8 . The 
luminosity bins have been chosen so that an equal 
number of objects are contained in each bin. The 
position of the data point along the luminosity axis 
indicates the average luminosity in this bin, while 
the error bars in luminosity indicate the range of 
luminosities covered. 

Based on the luminosity function, the contribu- 
tion of the AGNs to the total X-ray emissivity W 
can be estimated (Sazonov & Revnivtsev 2004). 

This can be done by simply multiplying the XLF 
by the luminosity in each bin and integrating 
over the range of luminosities ( 10 41 ergss _1 < 

-£' 20 — 40 keV < 10 45 5 ergs s" 1 ). This results in 
H / 20-4Okev(> 10 41 ergss" 1 ) = 2.8xl0 38 ergss ” 1 hf 0 Mpc" 3 . 
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5. Discussion 

The number flux distribution (Fig. 4) shows 
a slope of a = 1.66 ±0.11. In the local Uni- 
verse with no evolution, and isotropic distribution 
of AGN and assuming Eucledian geometry, the ex- 
pected value is a = 1.5 Even though the difference 
is only significant on a 1.5cr level, this might indi- 
cate that the area density at the low flux end of the 
distribution has been slightly over corrected. One 
has to keep in mind that only a few sources derived 
from a small area of the sky are constraining the 
low flux end. Krivonos et al. (2005) studied the 
extragalactic source counts as observed by INTE- 
GRAL in the 20-50 keV energy band in the Coma 
region. Based on 12 source detections they deter- 
mine a surface density of (1.4 ± 0.5) x 10 -2 deg -2 
above a threshold of 10“ 11 ergs cm" 2 s -1 in the 
20 — 50keV energy band, where we get a consis- 
tent value of (1.2 ± 0.2) x 10 -2 deg -2 . Comparing 
the total flux of all the objects in the AGN sam- 
ple (/20 — 40 keV = 2.6 x 10“ 10 ergs cm -2 s -1 sr -1 ) 
with the flux of the X-ray background as pre- 
sented by Gruber et al. (1999) shows that the 
INTEGRAL AGN account only for about 1% of 
the expected value. This is expected when tak- 
ing into account the high flux limit of our sample: 
La Franca et al. (2005) have shown that objects 
with /2-iofceV > 10 -n ergscm -2 s -1 contribute 
less than 1% to the CXB. This flux limit translates 
to the faintest flux in our sample of fio-wkeV — 
5.6 x 10 -12 ergs cm -2 s -1 for aT = 1.9 power law 
spectrum. 

We compared the emissivity per unit volume of 
our objects VF20-40 kev(> 10 41 ergss -1 ) = 2.8 x 
10 38 ergss -1 hf 0 Mpc -3 with that found in the 3- 
20 keV band. Assuming an average power law of 
r = 1.9, the extrapolated value is W 3 _ 2 okev(> 
10 41 ergss -1 ) = 6.7 x 10 38 ergs s -1 h 3 0 Mpc -3 , 

which is a factor of 1.3 larger than the value 
measured by RXTE (Sazonov & Revnivtsev 2004) 
and consistent within the la error. Extrapolating 
our result down to 2-10 keV shows an emissiv- 
ity which is lower by a factor of 1.2 compared 
to the one derived from the HEAO-l all-sky map 
(Miyaji et al. 1994) and also within the statistical 
uncertainties. 

The luminosity function derived from the IN- 
TEGRAL 20 - 40 keV AGN sample appears to be 
consistent with the XLF in the 2 — 20 keV range. 


A turnover in the XLF at 2 ^ 4.6 x 10 43 ergss -1 
is observed (Fig. 6). Below this luminosity also 
the fraction of absorbed AGNs starts to be larger 
than that of the unabsorbed ones, although the 
effect is significant only on a la level (Fig. 7). 
Both effects have been seen also in the 2 — lOkeV 
(Ueda et al. 2003; La Franca et al. 2005) and in 
the 3— 20keV band (Sazonov & Revnivtsev 2004). 
This implies that we do detect a similar source 
population as at lower energies. If a large pop- 
ulation of absorbed AGNs is dominating the cos- 
mic X-ray background at ~ 30 keV as indicated by 
HE AO 1 A-4 measurements (Gruber et al. 1999), 
and the source population is the same through- 
out the Universe those objects would have to 
have luminosities ±(20-40 keV) < 10 41 ergss -1 , as 
might be indicated by the larger fraction of ab- 
sorbed sources toward lower luminosities (Fig. 7). 
Even though it has to be taken into account that 
the low luminosity end of the XLF is based only 
on a small number of objects, below this lumi- 
nosity the distribution between active and nor- 
mal galaxies becomes blurred. Additionally ob- 
jects like Ultra-luminous X-ray sources (ULX) 
and star-forming galaxies could provide the nec- 
essary emission. One interesting case in this cat- 
egory is the detection of NGC 4395 with a lu- 
minosity of ±(20—40 keV) = I A x 10 40 ergs s -1 , 
consistent with measurements by XMM-Newton 
which showed ±(2-10 keV) = 1-5 x 10 40 ergss -1 
(Vaughan et al. 2005). The central engine of this 
galaxy has been classified as an “ultra-luminous” 
source, possibly associatied with an intermediate 
mass black hole with Mbh — (4*® x 1O 4 M 0 
(McHardy et al. 2005). In addition, NGC 4395 
harbors a ULX with ±( 2 -iokeV) = 10 39 ergss -1 at 
a distance of 2.9 kpc from the center of the galaxy 
(Colbert & Ptak 2002). 

Another scenario is that the population of X- 
ray emitting sources depends on redshift, i.e. that 
there is an evolution of population in time, indi- 
cating that the fraction of absorbed sources might 
be higher at larger redshifts, although it should 
be noted that the latter effect is not clearly de- 
tectable in the 2 - lOkeV range. The fraction 
of absorbed sources seems to depend on lumi- 
nosity (Ueda et al. 2003; Treister & Urry 2005), 
as seen also in the 20 — 40keV band (Fig. 8). 
But some studies come to the conclusion that 
there is no evolution of A r // (Ueda et al. 2004; 
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Treister & Urry 2005), while others find the frac- 
tion of absorbed sources increasing with redshift 
(La Franca et al. 2005). The latter also find that 
a combination of effects (the fraction of absorbed 
AGN decreases with the intrinsic X-ray luminos- 
ity, and increases with the redshift) can be ex- 
plained by a luminosity-dependent density evo- 
lution model. They further show that the lu- 
minosity function of AGNs like those presented 
here peaks at z ~ 0.7 while high luminosity 
AGN peak at z ~ 2. Unified models also pre- 
dict, depending on the applied model, a frac- 
tion of 0.6 — 0.7 for high flux low redshift AGN 
(Treister & Urry 2005). Worsley et al. (2005) 
examined Chandra and X MM- Newton deep fields 
and come to the conclusion that the missing CXB 
component is formed by highly obscured AGNs 
at redshifts ~ 0.5 — 1.5 with column densities of 
the order of fx = 10 23 — 10 24 cm' 2 . Evidence for 
this scenario is also found in a study of Chandra 
and Spitzer data (Polletta et al. 2006). Combin- 
ing multiwavelength data, this work estimates a 
surface density of 25 AGN deg -2 in the infrared 
in the 0.6 deg 2 Ch an draj SWIRE field, and only 
33% of them are detected in the X-rays down 
to /o. 3 - 8 keV = 10~ 15 ergs cm -2 s~U The work 
also indicates a higher abundance of luminous and 
Compton-thick AGN at higher redshifts (z > 0.5). 
This source population would be missed by the 
study presented here, because of the low redshifts 
(z = 0.022) of the INTEGRAL AGNs. 

Several studies (Ueda et al. 2004; Treister & Urry 
propose that the absorbed AGNs needed to ex- 
plain the CXB should be Compton thick, and 
therefore would have been missed at 2 — lOkeV. 
This argument does not hold for the INTEGRAL 
observations, where the impact of absorption is 
much less severe than at lower energies. The ef- 
fect on the measured flux of a source with pho- 
ton index T — 2 for Compton thick absorption 
(Nfj = 10 24 cm* 2 ) is only a 5% decrease in flux 
(40% for Nff = 10 25 cm -2 ). It is therefore un- 
likely that many Compton-thick objects have been 
missed by the INTEGRAL studies performed to 
date. One possibility would be, that they are 
among the newly detected sources found by IN- 
TEGRAL. The fraction of unidentified objects 
among those sources is of the order of 50%. It 
should be pointed out though, that most of these 
sources are located close to the Galactic plane and 


are more likely to belong to the Galaxy: the Sec- 
ond IBIS/ISGRI Soft Gamma-Ray Survey Cata- 
log (Bird et al. 2006) lists 55 new sources detected 
by INTEGRAL , of which 93% are located within 
-10° < b < +10°. Among those sources 3 are 
listed as extragalactic sources, 18 are of Galactic 
origin, and 29 have not been identified yet. 

In addition, those objects which have been 
demonstrated to be AGNs, are no more likely to 
be Compton-thick objects than the overall AGN 
population studied here. Only four AGNs (NGC 
1068, NGC 4945, MRK 3, Circinus galaxy) de- 
tected by INTEGRAL have been proven to be 
Compton thick objects so far, and none of them 
showed absorbtion of Nh > 5 x 10 24 cm~ 2 . In 
order to clarify this point, observations at soft X- 
rays of those objects without information about 
intrinsic absorption are required for all INTE- 
GRAL detected AGNs (Tab. 1). Up to now 23 
% of the INTEGRAL AGN are missing absorp- 
tion information. A first indication of what the 
absorption in these sources might be, can be de- 
rived from comparison of the INTEGRAL fluxes 
with ROSAT All-Sky Survey (RASS) Faint Source 
Catalogue data (Voges et al. 2000). In order to 
do so we assumed a simple power law with photon 
index F = 2.0 between the ROSAT 0.1 - 2.4keV 
band and the INTEGRAL 20 — 40keV range and 
fit the absorption. In the six cases where no detec- 
tion was achieved in the RASS, an upper limit of 
/(o. 1—2.4 keV) < 10 _13 ergscm“ 2 s“ 1 has been as- 
2005^umed, resulting in a lower limit for the absorp- 
tion Nh > (5 - 11) x 10 22 cm -2 . In Fig. 9 we show 
the distribution of intrinsic absorption. It has to 
be pointed out that the estimated values can only 
give an idea about the distribution of intrinsic ab- 
sorption and should not be taken literally, as the 
spectral slope between the measurements is un- 
known and the observations are not simultaneous. 
Nevertheless apparently none of the RASS detec- 
tions and non-detections requires an intrinsic ab- 
sorption of Nh > 2 x 10 23 cm -2 . Therefore it 
appears unlikely that a significant fraction of IN- 
TEGRAL AGNs will show an intrinsic absorption 
Nh > 10 24 cm* 2 . A similar result in the 3—20 keV 
band let Sazonov & Revnivtsev (2004) to the con- 
clusion that the missing emission to explain the 
CXB is not produced in ’normal’ AGNs, but that 
a comparable X-ray flux might be produced to- 
gether by low luminosity AGNs, non- active galax- 
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ies and clusters of galaxies. 

Most investigations have been focused so far 
on the X-rays below 20 keV, and INTEGRAL 
can add substantial information to the nature of 
bright AGNs in the local Universe. Considering 
the expected composition of the hard X-ray back- 
ground, it does not currently appear that the pop- 
ulation detected by INTEGRAL can explain the 
peak at 30 keV, as Compton thick AGNs are ap- 
parently less abundant than expected. But the 
sample presented here might be still too small to 
constrain the fraction of obscured sources, and 
the missing Compton thick AGN could be de- 
tected when studying sources with /(20-40keV) < 
10“ n ergscm _2 s“ 1 . 

6. Conclusions 

A statistically complete extragalactic sample 
derived from the INTEGRAL public data archive 
comprises 58 low redshift Seyfert galaxies (( z ) = 
0.022 ± 0.003) and 8 blazars in the hard X-ray 
domain. Two galaxy clusters are also detected, 
but no star-burst galaxy has been as yet. This 
INTEGRAL AGN sample is thus the largest one 
presented so far. 

The number flux distribution shows a slope 
of a = 1.66 ±0.11. Because of the high flux 
limit of our sample the objects account in to- 
tal for less than 1% of the 20 — 40keV cosmic 
X-ray background. The emissivity of all AGNs 
per unit, volume U / 2 o- 4 okev(> 10 41 ergs s" 1 ) = 
2.8 x 10 38 ergss -1 h 3 0 Mpc” 3 appears to be con- 
sistent with the background measurements in the 
2-10 keV energy band. 

The luminosity function in the 20 — 40 keV en- 
ergy range is consistent with that measured in the 
2 - 20keV band. Below the turnover luminosity 
of L* = 4.6 x 10 43 ergss _1 the absorbed AGNs be- 
come dominant over the unabsorbed ones. Only 
a small fraction of Compton thick AGNs are de- 
tectable within the sample with known intrinsic 
absorption. For the sources without reliable ab- 
sorption information we derived an estimate from 
the comparison with ROSAT All-Sky Survey data 
and find that the data do not require additional 
Compton thick objects within the sample pre- 
sented here. The sources to explain the 30 keV 
peak in the cosmic X-ray background have to be 
either at higher redshift, as proposed for example 


by Worsley et al. (2005), which would imply an 
evolution towards stronger absorption with red- 
shift, or the missing emission has to be due to low 
luminosity AGNs, clusters of galaxies, non-active 
galaxies or ultra-luminous X-ray sources. 

Over the life time of the INTEGRAL mission 
we expect to detect of the order of 200 AGNs. 
Combining these data with the studies based on 
Swift /BAT, operating in a similar energy band as 
IBIS/ISGRI, will further constrain the hard X-ray 
luminosity function of AGNs. But we will still be 
limited to relatively the high flux, low luminosity 
and low redshift end of the distribution, which will 
be inadequate to explain the cosmic X-ray back- 
ground at E > 20keV. Future missions such as 
EXIST or NuStar will be required to answer the 
question of what dominates the Universe in the 
hard X-rays. 
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Fig. 1. — Exposure map representing the data 
used in our analysis in Galactic coordinates. Con- 
tours indicate 2 ks, 5 ks, 10 ks, 20 ks, 100 ks, 
200 ks, 500 ks, 1 Ms, 2 Ms exposure time. INTE- 
GRAL spent most of the observing time on and 
near the Galactic plane, with a strong focus on the 
Galactic center and on areas including bright hard 
X-ray sources like the Cygnus region, Vela, GRS 
1915+105, and the Crab. Fields at high Galac- 
tic latitude include Coma, Virgo and the region 
around NGC 4151. 



log t [sec] 

Fig. 2. — Correlation of exposure time and flux 
limit (5a, 20 — 40keV) for the AGNs in this study. 
The curve shows a smooth polynomial fit for flux 
limit versus logarithm of the exposure time. 
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Table 1 

INTEGRAL AGN catalog 


Name 

Type 

z 

R.A. 

(J2000.0) 

Decl. 

(J2000.0) 

exp. a 

(ks) 

ISGRI 

f x b 

log 7/20-40 keV 

(ergs s _1 ) 

N h cJ 

Ref. f 

1ES 0033+595 

BL Lac 

0.086 

00 35 53 

+59 50 05 

1449 

3.5 

0.37 

43.83 

0.36 1 

2 

NGC 788 

Sy 1/2 

0.0136 

02 01 06 

-06 48 57 

311 

10.7 

2.98 

43.09 

< 0.02 3 

4 

NGC 1068 

Sy 2 

0.003793 

02 42 41 

-00 00 48 

311 

4.3 

0.93 

41.47 

> 150 s 

4 

QSO B0241+62 

Sy 1 

0.044557 

02 44 58 

+62 28 07 

43 

3.4 

2.02 

43.97 

1.5 6 

7 

NGC 1142 

Sy 2 

0.028847 

02 55 12 

-00 11 02 

311 

5.5 

1.58 

43.48 

? 

8 

NGC 1275 

Sy 2 

0.017559 

03 19 48 

+41 30 42 

264 

8.4 

1.89 

43.12 

3.75 4 

4 

3C 111 

Sy 1 

0.048500 

04 18 21 

+38 01 36 

67 

10.0 

6.27 

44.54 

0.63 3 

4 

UGC 3142 

Sy l 

0.021655 

04 43 47 

+28 58 19 

247 

16.8 

5.46 

43.76 

*? 

2 

LEDA 168563 

Sy 1 

0.0290 

04 52 05 

+49 32 45 

28 

2.8 

2.27 

43.64 

? 

2 

MCG +8-11-11 

Sy 1.5 

0.020484 

05 54 54 

+46 26 22 

21 

6.2 

6.07 

43.76 

< 0.02 6 

4 

MRK 3 

Sy 2 

0.013509 

06 15 36 

+74 02 15 

472 

15.9 

3.65 

43.17 

no 6 

4 

MRK 6 

Sy 1.5 

0.018813 

06 52 12 

+74 25 37 

482 

8.7 

2.01 

43.21 

10 6 

4 

S5 0716+714 

BL Lac 

0.3 d 

07 21 53 

+71 20 36 

482 

0.7 

0.14 

44.41 d 

< 0.01 3 

9 

ESO 209-12 

Sy 1.5 

0.040495 

08 01 58 

-49 46 36 

1543 

6.7 

0.86 

43.52 

? 

7 

FRL 1146 

Sy 1 

0.031789 

08 38 31 

-35 59 35 

849 

3.6 

0.60 

43.15 

? 

7 

S5 0836+710 

FSRQ 

2.172 

08 41 24 

+70 53 42 

391 

6.4 

1.73 

47.79 

0.11 3 

9 

MCG-05-23-1 6 

Syl.9 

0.008486 

09 47 40 

-30 56 56 

2 

2.3 

11.20 

43.25 

1.6 16 

16 

IGR J 10404-4625 

Sy 2 

0.0237 

10 40 22 

-46 25 26 

46 

1.5 

0.67 

42.93 

7 

10 

IGR J 12026-5349 

AGN 

0.028 

12 02 48 

-53 50 08 

191 

5.5 

1.86 

43.52 

2.2 11 

11 

NGC 4051 

Sy 1.5 

0.002336 

12 03 10 

+44 31 53 

443 

8.4 

1.80 

41.34 

< 0.01 6 

4 

NGC 4151 

Sy 1.5 

0.003320 

12 10 33 

+39 24 21 

483 

163.3 

26.13 

42.80 

6.9 12 

12 

NGC 4253 

Sy 1.5 

0.012929 

12 18 27 

+29 48 46 

715 

6.1 

0.93 

42.54 

0.8 6 

4 

4C +04.42 

BL Lac 

0.965 

12 22 23 

+04 13 16 

690 

4.5 

0.80 

46.58 

? 

7 

NGC 4388 

Sy 2 

0.008419 

12 25 47 

+ 12 39 44 

215 

34.8 

9.54 

43.18 

27 13 

13 

NGC 4395 

Sy 1.8 

0.001064 

12 25 49 

+33 32 48 

739 

5.1 

0.56 

40.14 

0.15 3 

4 

3C 273 

Blazar 

0.15834 

12 29 07 

+02 03 09 

655 

34.2 

5.50 

45.58 

0.5 4 

14 

NGC 4507 

Sy 2 

0.011801 

12 35 37 

-39 54 33 

152 

14.9 

6.46 

43.30 

29 6 

4 

IGR J12391-1612 

Sy 2 

0.036 

12 39 06 

-16 10 47 

83 

1.4 

3.46 

44.02 

1.9 11 

11 

NGC 4593 

Sy 1 

0.009000 

12 39 39 

-05 20 39 

723 

20.1 

3.31 

42.78 

0.02 6 

4 

IGR 312415-5750 

Sy 2 

0.024 

12 41 24 

-57 50 24 

440 

1.1 

0.33 

42.64 

? 

15 

3C 279 

Blazar 

0.53620 

12 56 11 

+05 47 22 

497 

3.6 

0.82 

45.97 

< 0.13 3 

4 

Coma Cluster 

GClstr 

0.023100 

12 59 48 

+27 58 48 

516 

7.2 

1.09 

43.11 

< 0.01 4 

4 

NGC 4945 

Sy 2 

0.001878 

13 05 27 

-49 28 06 

276 

33.8 

9.85 

41.88 

400 6 

16 

ESO 323-G077 

Sy 2 

0.015014 

13 06 26 

-40 24 53 

761 

6.9 

1.20 

42.78 

55 17 

17 

IGR J13091 + 1137 

AGN 

0.0251 

13 09 06 

+ 11 38 03 

48 

2.0 

1.06 

43.18 

90 11 

15 

NGC 5033 

Sy 1.9 

0.002919 

13 13 28 

+36 35 38 

377 

4.6 

1.06 

41.30 

2.9 6 

4 

Cen A 

Sy 2 

0.001830 

13 25 28 

-43 01 09 

532 

167.4 

32.28 

42.38 

12. 5 4 

16 

MCG-06-30-015 

Sy 1 

0.007749 

13 35 54 

-34 17 43 

567 

4.9 

0.73 

41.99 

7.7 6 

4 

4U 1344-60 

AGN 

0.043 

13 47 25 

-60 38 36 

603 

16.6 

2.83 

44.09 

2.19 6 

4 

IC 4329A 

Sy 1.2 

0.016054 

13 49 19 

-30 18 36 

440 

41.7 

8.19 

43.68 

0.42 6 

4 

Circinus gal. 

Sy 2 

0.001448 

14 13 10 

-65 20 21 

589 

58.9 

10.73 

41.69 

360 6 

16 

NGC 5506 

Sy 1.9 

0.006181 

14 13 15 

-03 12 27 

101 

12.6 

4.21 

42.55 

3.4 6 

4 

NGC 5548 

Sy 1.5 

0.017175 

14 18 00 

+25 08 12 

211* 

2.4 

0.71 

42.67 

0.51 6 

4 

PG 1416-129 

Sy l 

0.129280 

14 19 04 

-13 10 44 

117 

8.3 

4.86 

45.33 

0.09 4 

4 

ESO 511-G030 

Sy l 

0.022389 

14 19 22 

-26 38 41 

145 

0.1 

1.93 

43.34 

C0.05 17 

17 

IC 4518 

Sy 2 

0.015728 

14 57 43 

-43 07 54 

338 

2.3 

0.49 

42.44 

? 

4 

IGR 316119-6036 

Sy l 

0.016 

16 11 54 

-60 36 00 

475 

1.2 

0.25 

42.16 

? 

2 

IGR J1 6482-3036 

Sy 1 

0.0313 

16 48 17 

-30 35 08 

973 

4.2 

0.73 

43.22 

? 

10 

NGC 6221 

Sy 1/2 

0.004977 

16 52 46 

-59 13 07' 

523 

5.6 

1.32 

41.86 

I 18 

4 

Oph Cluster 

GClstr 

0.028 

17 12 26 

-23 22 33 

1763 

30.8 

4.10 

43.90 

? 

25 

NGC 6300 

Sy 2 

0.003699 

17 17 00 

-62 49 14 

173 

10.0 

3.91 

42.07 

22 19 

16 

GRS 1734-292 

Sy 1 

0.021400 

17 37 24 

-29 10 48 

3332 

45.9 

4.03 

43.62 

3.7 4 

24 

2E 1739.1-1210 

Sy 1 

0.037 

17 41 54 

-12 11 52 

631 

5.5 

1.03 

43.51 

? 

7 

IGR J 1 8027-1455 

Sy 1 

0.035000 

18 02 47 

-14 54 55 

942 

12.6 

2.03 

43.76 

19. 0 4 

20 

PKS 1830-211 

Blazar 

2.507 

18 33 40 

-21 03 40 

1069 

12.7 

2.07 

48.02 

< 0.7 3 

21 

ESO 103-G35 

Sy 2 

0.013286 

18 38 20 

-65 25 39 

36 

4.2 

2.97 

43.07 

15 3 

16 

2E 1853.7+1534 

Sy 1 

0.084 

18 56 00 

+ 15 38 13 

761 

9.3 

1.74 

44.48 

? 

10 
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Fig. 3. — The total survey area with respect 
to the 5 a flux limit in the 20 - 40 keV band. 
For a flux limit of fx > 10" 11 ergs cm” 2 s~ l the 
survey covers 76% of the sky, and for fx > 
10 _11 ergscm _2 s _1 17%. 



- 1.5 - 1.0 - 0.5 0.0 0.5 1.0 

S (Flux Threshold, x 1 CT’Vgs/em 1 /*) 

Fig. 4. — Number flux distribution per steradian 
of INTEGRAL AGN with a detection significance 
> 5cj. Blazars have been excluded. The maxi- 
mum likelihood slope as described in Section 3.1 
is 1.66 ±0. 11- 



Fig. 5. — Cumulative luminosity function of IN- 
TEGRAL AGN with a detection significance > 
5cr. Blazars have been excluded. The density <t> 
describes the number of objects per Mpc 3 above a 
given luminosity Lx • 
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Fig. 6. — Differential luminosity function of AGN 
with A log Lx = 0.5 binning. The line shows a fit 
to a smoothly connected two power-law function 
with a turnover luminosity at 4.6 x 10 43 ergss” 1 . 
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Fig. 7. — Cumulative AGN luminosity function 
for 21 absorbed ( Nh > 10 22 cm~ 2 ; triangles) and 
13 unabsorbed sources (octagons). 
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Fig. 9. — Distribution of intrinsic absorption for 
all INTEGRAL AGN (blazars excluded), as mea- 
sured in the soft X-rays. The shaded area shows 
the reliable measurements, the other values are 
based on comparison of ROSAT All-Sky Survey 
and INTEGRAL data. Lower limits on absorp- 
tion have been excluded. 


Fig. 8. — Fraction of absorbed AGNs (Nh > 
10 22 cm" 2 ) versus luminosity. The position of the 
data point along the luminosity axis indicates the 
average luminosity in this bin, while the error bars 
in luminosity indicate the range of luminosities 
covered. 
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Table 1 — Continued 


Name 

Type 

z 

R.A. 

(J2000.0) 

Decl. 

(J2000.0) 

exp. a 

(ks) 

ISGRI 

( CT ) 

fx b 

log L 20-40 keV 
(ergs s" 1 ) 

N„ c -' 

Ref. 1 

1H 1934-063 

Sy 1 

0.010587 

19 37 33 

-06 13 05 

684 

2.7 

0.48 

42.08 

7 

4 

NGC 6814 

Sy 1.5 

0.005214 

19 42 41 

-10 19 25 

488 

12.1 

2.92 

42.24 

< 0.05 3 

4 

IGR J19473+4452 

Sy 2 

0.0539 

19 47 19 

+44 49 42 

969 

5.9 

1.05 

43.86 

ll 11 

11 

Cygnus A 

Sy 2 

0.056075 

19 59 28 

+40 44 02 

1376 

21.6 

3.24 

44.39 

20 22 

4 

MCG +04-48-002 

AGN 

0.014206 

20 28 35 

+25 44 00 

187 

3.1 

1.10 

42.70 

7 

2 

4C +74.26 

AGN 

0.0104 

20 42 37 

+75 08 02 

72 c 

1.9 

0.93 

42.35 

0.21 3 

23 

MRK 509 

Sy 1 

0.034397 

20 44 10 

-10 43 25 

73 

8.6 

4.66 

44.11 

< 0 . 01 3 

4 

IGR J21247+5058 

AGN 

0.020 

21 24 39 

+50 58 26 

213 

11.9 

4.15 

43.57 

7 

20 

3C 454.3 

Blazar 

0.859 

22 53 58 

+16 08 54 

92“ 

2.8 

3.1 

46.56 

7 

26 

MR 2251-178 

Sy i 

0.063980 

22 54 06 

-17 34 55 

489 

7.0 

1.20 

44.07 

< 0.19 3 

4 

MCG -02-58-022 

Sy 15 

0.046860 

23 04 44 

-08 41 09 

489 

3.9 

1.20 

43.79 

< 0.08 3 

4 


a ISGRI exposure time 
b flux /(20— 40 keV) in 10 -n ergcm^s" 1 
“intrinsic absorption in 10 22 cm -2 
tentative redshift 

“not covered by survey presented here 

f REFERENCES. — (1) Donato. Sambruna, Gliozzi 2005; (2) Bassani et al. 2006; (3) Tartarus database; (4) Beckmann et a!. 2006a; 
(5) Matt et al. 1997; (6) Lutz et al. 2004; (7) Bird et al. 2006; (8) Virani et al. 2005; (9) Pian et al. 2005; (10) Masetti et al. 2005; 
(11) Sazonov et al. 2005; (12) Beckmann et al. 2005; (13) Beckmann et al. 2004; (14) Courvoisier et al. 2003b; (15) Revnivtsev et 
al. 2005; (16) Soldi et al. 2005; (17) Sazonov & Revnivtsev 2004; (18) Levenson, Weaver, fe Heckman 2001; (19) Matsumoto, Nava, 
Maddox et al. 2004; (20) Masetti et al. 2004; (21) De Rosa et al. 2005; (22) Young et al. 2002; (23) This work; (24) Sazonov et al. 
2004; (25) Revnivtsev et al. 2004; (26) Pian et al. 2006 


14 


